191 Search Results

Continuum computational kinetic plasma models evolve the distribution function of a plasma species f_s on a phase-space grid over time. In many problems of interest the distribution function has limited extent in velocity space; hence, using a uniform, highly refined mesh would be costly and slow. Nonuniform velocity grids can reduce the computational cost by placing more degrees of freedom where f_s is appreciable and fewer where it is not. In this work we introduce a first-of-its kind discontinuous Galerkin approach to nonuniform velocity-space discretization using mapped velocity coordinates. This new method is presented in the context of a gyrokineticmore » model used to study magnetized plasmas. We create discretizations of collisionless and collisional terms using mappings in a way that exactly conserves particles and energy. Numerical tests of such properties are presented, and we show that this new discretization can reproduce earlier gyrokinetic simulations using grids with up to 6–60 times fewer cells and 22X-60X speed-ups depending on dimensionality, geometry and plasma parameters.« less

Land surface relative humidity (RH) is a key variable in the coupled land–atmosphere system that profoundly influences terrestrial hydroclimate and ecosystems. Yet historical changes in land RH are not well understood due to limited observations, biased reanalyses, and the lack of a framework for interpreting RH changes under multiple influencing factors. Here, we show that the spatiotemporal variability of land RH and its distinct historical trends among observations, reanalyses, and Earth system models are captured by a simple index based on the ratio of precipitation (P) to a modified potential evapotranspiration formulated independently of RH (PET°). The index provides amore » physical calibration of biased land RH in reanalyses and a quantitative framework for interpreting land RH changes. Over 1973–2024, land RH has decreased substantially, owing to the intrinsic rise in PET° with temperature and little increase in land precipitation. Reanalyses overestimate the observed RH decrease, consistent with exaggerated surface warming and precipitation decline. The index captures this coherent bias and enables a calibration using observed precipitation and temperature. Models simulate a wide range of land RH trends, but nearly all runs underrepresent the historical drying. The index captures the model spread and discrepancy and attributes them to contributions of precipitation and PET°. Weaker land RH decreases in models arise mainly from weaker subtropical precipitation declines, linked to muted intensification of subtropical highs and biased subtropical climatology. The model–observation discrepancy is unlikely explained by internal variability, implying model underestimation of forced RH decrease and a drier land future than current projections.« less

We report the in situ formation of Ni nitride for plasma-assisted ammonia synthesis. Both the surface nitrogen concentration and the ammonia formation rate exhibit dependence on the N₂:H₂ feed ratio. The maximum surface nitrogen concentration occurs at a N₂:H₂ ratio of 4:1, and the maximum catalytic activity occurs at 2:1. In contrast, the formation of gas phase radicals is less sensitive to feed composition, indicating that Ni nitride is more kinetically relevant to ammonia production than gas-phase radicals. The plasma-induced formation of Ni nitride is therefore proposed to be a critical contributor to the synergistic effects in plasma-assisted catalytic ammoniamore » synthesis. Additionally, Ni nitride alters the surface reaction mechanism of plasma-assisted ammonia synthesis, with the rate-determining-step (RDS) shifting to surface-bound NH₃ formation rather than N₂ activation at temperatures below 373 K. These findings provide mechanistic insight that opens opportunities for optimizing the performance of plasma-assisted catalytic ammonia synthesis« less

Recent separated reactant experiments for thin-shell (6 µ⁢m) shock-driven implosions on OMEGA have demonstrated significant mix from a buried deuterated layer of the shell into the hot spot. Time resolved D³He-p reaction history data demonstrate a (50 ± 20)⁢ ps shift earlier in peak nuclear emission for separated reactant experiments relative to control, in contrast to past experimental data for thicker, 20 µ⁢m shells with no laser burn through that show a 75 ps delay due to the time required for hydrodynamic instabilities to develop. This contrast suggests that the mix mechanism was not hydrodynamic. Ion kinetic simulations utilizing fallmore » line analyses show much closer agreement with mix yield and temperature than diffusion models, predicting a D³He-p mix yield of 1.7 × 10⁹ as compared to the experimental value of 9.3⁢ (±2.1) × 10⁸. This is three orders of magnitude closer than the fall line analysis from a hydrodynamic simulation with an inline diffusive mix model, which suggests minimal mix and D³He-p yields of 5×10⁵. This makes kinetic mechanisms the only feasible explanation for the mix seen, demonstrating impact of a non-standard mix mechanism. An analytical model of this kinetic mix mechanism suggests that it can remain significant in situations when the shell expands significantly to low densities, and diffusive models predict negligible mix. Finally, kinetic mix will impact multiple types of high energy density, laser-driven fusion experiments including high-adiabat direct drive cryoexperiments, nuclear cross section experiments, and thin-shell polar direct drive experiments used to tune heat conduction models.« less

The edge turbulence model Hermes (Dudson et al 2017 Plasma Phys. Control. Fusion 59 05401) is set up for plasma boundary simulations with an radiofrequency (RF) antenna, using parameters characteristic of a tokamak edge. Cartesian slab geometry is used with thin plate limiters representing the ion cyclotron range of frequency (ICRF) antenna side-wall limiters. Ad-hoc DC electric biasing of the limiters, motivated by calculations with VSim (Nieter et al 2004 J. Comput. Phys. 196 448), represents an induced RF sheath rectified potential in the plasma turbulence model. Flux-driven turbulence simulations demonstrate a realistic distribution of plasma profiles and fluctuations. Theremore » is a clear effect of the antenna sheath voltage leading to formation of convective cells; bias-induced convective transport flattens the scrape-off layer density profile and fluctuations penetrate into the shadow region of the limiters as the bias voltage increases. Turbulent transport for impurity ions is inferred by following ion trajectories in the simulated plasma turbulence fields, showing Bohm-like effective diffusion rates. All in all, the model elucidates the key physical phenomena governing the effects of ICRF-induced antenna biasing on the tokamak boundary plasma.« less

Removing organic ligands from colloidal nanoparticles is critical for fabricating solid-state devices, yet accurately quantifying this removal remains a significant analytical challenge. Here, we establish a robust and accessible method for this quantification by calibrating Raman spectroscopy against precise ion beam analysis (IBA) for nanoparticle assemblies (CNAs) processed by helium (He) and oxygen (O₂) plasmas. We demonstrate that the calibration curves are remarkably independent of plasma power and pressure, depending critically only on the choice of feed gas. He plasma induces rapid dehydrogenation and cross-linking, evidenced by a much faster decrease in the C–H Raman signal relative to the actualmore » carbon loss. Conversely, O₂ plasma leads to a surprising “apparent hydrogenation”, where the carbon backbone is removed significantly faster than the C–H signal diminishes. This counterintuitive effect is explained by a serial mechanism of oxidative fragmentation; β-scission cleaves the alkyl chains, and subsequent stabilization steps enrich the remaining film with hydrogen-rich methyl-terminated fragments, while carbon is efficiently removed as volatile CO. This work provides calibrated functions that enable the rapid determination of absolute carbon content in processed CNAs using simple Raman spectroscopy with uncertainties of ∼8% for O₂ and ∼12% for He plasma, offering a vital tool for both process diagnostics and fundamental studies of plasma–matter interactions in colloidal nanocrystal films.« less

Ru₃Sn₇ is experimentally demonstrated as an efficient catalyst, with potential utilization of topological surface states for hydrogen evolution reaction. Despite its promising catalytic performance, the topological nature of Ru₃Sn₇ remains uncertain. Particularly, the bulk-boundary correspondence has not yet been established, hence hindering a rigorous justification of its topologically-protected surface states. In this work, the bulk topology of Ru₃Sn₇ is detailed using first-principles calculations and the topological quantum chemistry formalism. Ru3Sn7 turns out to be an enforced semimetal possessing symmetry-protected crossings within a set of bands near the Fermi level, which are enforced and prescribed by the violations of symmetry-prescribed compatibilitymore » relations. Moreover, the surface states and the associated origin from the same set of entangled bands are identified, thereby establishing the bulk-boundary correspondence. To evaluate the effects of chemical modifications, the response of topological surface states to various surface terminations, stoichiometry, and oxidation is examined. The surface structures are globally optimized, and the phase diagrams for various experimental conditions are built. It is shown that, due to changes in the local chemical environment, the original surface states are significantly altered. Modified surface bands can be observed near the Fermi level on surface terminations that preserve the C_4v symmetry.« less

The scrape-off layer of a tokamak fusion reactor carries the plasma exhaust from the hot core plasma to the material surfaces of the reactor vessel. The heat loads imposed by the exhaust are a critical limit on the performance of fusion power plants. Turbulent transport of the plasma regulates the width of the scrape-off layer plasma and must be modelled to understand the intensity of these heat loads. STORM is a plasma turbulence code capable of simulating three dimensional turbulence across the full scrape-off layer of a tokamak fusion reactor, using a drift reduced, collisional fluid model. STORM uses mostlymore » finite difference schemes, with a staggered grid in the direction parallel to the magnetic field. We describe the model, geometry and initialisation options used by STORM, as well as the numerical methods, which are implemented using the BOUT++ plasma simulation framework. BOUT++ has been enhanced alongside the development of STORM, providing better support for staggered grid methods. We summarise these enhancements, including a detailed explanation of the parallel derivative methods, which underwent a major update for version 4 of BOUT++.« less

A velocity dependent, kinetic model for line radiation is developed for continuum kinetic codes. It has been implemented in the full-f gyrokinetic code Gkeyll. The total radiation for a charge state is modeled as an advection in velocity space with a form of $$\nabla_v \cdot(v\nu(v)f(v))$$, guaranteeing particle conservation. The velocity dependence (in the form of an effective frequency $$\nu(v)$$) is found through fitting the energy loss of the operator, i.e. the second velocity moment, to the radiation data in the OpenADAS database. Therefore, each individual transition does not need to be evaluated every time step, significantly reducing the computational costmore » of including line radiation in a kinetic model. The dependence on velocity instead of the usual, temperature, allows the radiation to be computed from non-Maxwellian electron distribution functions: We benchmark the model against a collisional radiative model using isotropic non-Maxwellian distribution functions. A velocity dependent model of radiation can more accurately describe the radiation in the more kinetic regimes expected in reactor-scale devices. The velocity dependence qualitatively captures the quantum mechanical need for a minimum velocity before any radiation occurs.« less

This paper presents the first gyrokinetic (GK) simulations of edge and scrape-off layer (SOL) turbulence accelerated by a velocity-space spectral approach in the full-f GK code GENE-X. Building upon the original grid velocity-space discretization, we derive and implement a new spectral formulation and verify the numerical implementation using the method of manufactured solution. We conduct a series of spectral turbulence simulations focusing on the TCV-X21 reference case (Oliveira et al., 2022 [26]) and compare these results with previously validated grid simulations (Ulbl et al., 2023 [25]). The spectral approach reproduces the outboard midplane (OMP) profiles (density, temperature, and radial electricmore » field), dominated by trapped electron mode (TEM) turbulence, with excellent agreement and significantly lower velocity-space resolution. As a consequence, the spectral approach reduces the computational cost (CPUh) by at least an order of magnitude, of approximately 50 for the TCV-X21 case. This enables high-fidelity GK simulations to be performed within a few days on modern CPU-based supercomputers for medium-sized devices and establishes GENE-X as a powerful tool for studying edge and SOL turbulence, moving towards reactor-relevant devices like ITER.« less